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SOME PROPERTIES OF EVEN MOMENTS OF 
UNIFORM RANDOM WALKS 


KEVIN G. HARE AND GHISLAIN MCKAY 


Abstract. We build upon previous work on the densities of uni¬ 
form random walks in higher dimensions, exploring some proper¬ 
ties of the even moments of these densities and extending a result 
about their modularity. 


1. Introduction 


Consider a short random walk of n steps in d dimensions where 
each step is of unit length and whose direction is chosen uniformly. 
Following [2], we let 1 / = I — 1 and denote by pn{i^',x) the probability 
density function of the distance x to the origin of this random walk. 
This paper will be concerned with the even moments of these random 
walks. 

Definition 1.1. Define 



as the moment of the probability density function. 

We know that 

Theorem 1.2 (Borwein, Staub, Vignot, Theorem 2.18, [ 2 ]). For non¬ 
negative integers k, Wn{m, 2k) is given by 



k nu 

ki + v, ... ,kn + v 


fclH-hfcn —k 


Theorem 1.3 (Borwein, Staub, Vignot Example 2.23, [2]). For given 
V, let A{i>) be the infinite lower triangular matrix with entries 
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for row indices k = 0,1,2,... and columns entries j = 0,1,2,.... 
Then the moments Wn+i{i''-,‘2^k) are given by the row sums of . 


For a good history of these moments, and random walks in general, 

see [HEIEIII]. 


Example 1.4. For example, the upper corner of A{0),A{1) and A{2) 
are given below. 


A(0) := 




A(2) := 


1 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

1 

4 

1 

0 

0 

0 

0 

0 

1 

9 

9 

1 

0 

0 

0 

0 

1 

16 

36 

16 

1 

0 

0 

0 

1 

25 

100 

100 

25 

1 

0 

0 

1 

36 

225 

400 

225 

36 

1 

0 

1 

49 

441 

1225 

1225 

441 

49 

1 


■ 1 

0 

0 

0 

0 

0 

0 0 

1 

1 

0 

0 

0 

0 

0 0 

1 

3 

1 

0 

0 

0 

0 0 

1 

6 

6 

1 

0 

0 

0 0 

1 

10 

20 

10 

1 

0 

0 0 

1 

15 

50 

50 

15 

1 

0 0 

1 

21 

105 

175 

105 

21 

1 0 

1 

28 

196 

490 

490 

196 

28 1 

1 

0 

0 

0 

0 

0 

0 0 

1 

1 

0 

0 

0 

0 

0 0 

1 

8/3 

1 

0 

0 

0 

0 0 

1 

5 

5 

1 

0 

0 

0 0 

1 

8 

15 

8 

1 

0 

0 0 

1 

35/3 

35 

35 

35/3 

1 

0 0 

1 

16 

70 

112 

70 

16 

1 0 

1 

21 

126 

294 

294 

126 

21 1 


The lower triangular entries of A(0) are the squares of the binomial 
coefficients (^.) and those in A(l) are known as the Naryana numbers 
m A001263]. Using these observations about A(0) and A(l), it is easy 
to observe that all of the coefficients of A(0) and A(l) are integers. A 
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quick glance at ^(2) shows that this is not always true. It was stated 
that Akj{2) G in [2]. 

We dehne 



Using this notation we see that ro = ri = 1 and r 2 = 3. It is not 
immediately that is well dehned and hnite for all u, (although we 
will show that this is the case). 

In Section [2] we show that 

Theorem 1.5. For u > 1 we have | 

This is not best possible. In Sectionj^we prove the opposite direction 

Theorem 1.6. For u > 1 we have I 

We conjecture that this is in fact best possible. That is, we conjecture 

Conjecture 1.7. For u > 1 we have r^ = 

We present evidence for this conjecture in Section and 
Next we consider a result by Borwein, Nuyens, Straub and Wan in 
[1] about the modularity of moments. They showed that 

Theorem 1.8. For primes p, we have 


Wn{0; 2p) = n mod p. 


We extend this in Section to get 

Theorem 1.9. Let 

• p = k be prime with 2z/ < p, or 

• p = k + u be prime with v < p. 

Then 


Wn{i^; 2k) = n mod p. 


If p‘^ = k with p prime then 


IW(0; 2k) = n mod p^. 


It is worth remarking that if both pi ■= k and p 2 '.= k + u are prime 
with 2z/ < Pi (and hence p < 2p < pi < P2), then clearly Wn{p] 2k) = n 
mod P 1 P 2 by the Chinese Remainder Theorem. 

In Section we discuss some of the open problems related to this 
research. 
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2. A PROOE OE Theorem 1.5: r^,\{2v - ly./ul 


To prove Theorem |1.5[ we make use of the following remark and 
lemma: 


Remark 2.1. There are multiple equivalent ways of representing . 

The three most common that we will use are: 

. . . _ fk\ {k + z/)!z/! 

-QiTT.T 


Lemma 2.2. For integers 1 < p < j we have 

gcd((j - z/ + l)(j - 1 / + 2) ■ ■ ■ j, (j + l)(j + 2) ■ ■ ■ (j + z/)) I (2z/ - 1)! 

Proof. Let Aj^i, = j — u + 1,... ,j and = j + 1,..., j + u. Let 
7r(Aj^i,) and be the products of these sequences. Let p be a 

prime number and Vp{x) be the p-adic valuation of x. We see that for 
p" > 2z/ that there is at most one term in Aj^y U Bj^y that is divisible 
by p". Without loss of generality we may assume that such a term, if 
it exists, is in Aj^^. We see that Vp{Bj^^) = Vp{Bj^pai. ^f) for all k by 
translation. Further, if there exists a term in Aj^y that is divisible by 
p", then, by translations we can assume that this term is divisible by 
an arbitrarily high power of p. Hence we can assume that, if such a 
term exists, then we can hnd a translate of this sequence so that 


5 ^(-®j+P“Zc,^))) jj^pak^uf) • 

We see that ii p^ < v then there are at most ^ terms in B 
are are divisible by p^. We see that if z/ < p^ < 2z/ then there are at 


j-\-p^k^u 


most 


2v 

HP 


1 terms in Bjj^pc,k,v are are divisible by p^. By Chinese 

remainder theorem we can hnd such a j so that both the inequalities 
are exact. This gives us that 

'2z/' 


(1) np(gcd(7r(Aj_|_pQ!fc J/), 7r(Rj_|_pa^ ^ 'y ^ 

p^<u 


V 

p/3 


E 

v<p^<2v 


pP 


-1 


and moreover there exists a j so that this is exact. 
We observe that 

Vp{{2u-iy.)= 

p/3<2!y-l 
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1 / 

Equation (|T|) 

(21.-1)! 

1 

1 

1 

2 

2-3 

2 ■ 3 

3 

23.3.5 

23.3 ■ 5 

4 

23.32 ■ 5 ■ 7 

2W 32.5.7 

5 

2® ■ 33 ■ 5 ■ 7 

2U 34.5.7 

6 

2® ■ 33 ■ 52 ■ 7 ■ 11 

2® ■ 3^ ■ 52 ■ 7 ■ 11 

7 

2U 34.52.7.11. 13 

210.35.52.7.11. 13 

8 

27.34.52.72.11. 13 

211.3® ■ 53 ■ 72 ■ 11 ■ 13 

9 

211 . 34.52.72.11 . 13 . 17 

215.3® ■ 53 ■ 72 ■ 11 ■ 13 ■ 17 

10 

241 . 36.52.72.11 . 13.17.19 

216.3® ■ 53 ■ 72 ■ 11 • 13 ■ 17 ■ 19 


Table 1. Prime factorization of Eq Q and (2i/ — 1)!. 


Observe that ii < u then 

2u-l 


P‘ 


i/3 


> 


P' 


If p/3 = j/ then 


li 1 / < < 21 / — 1 then 


2i/-l 


1 / 

1 

_1 


pP 


= 1 . 


Lastly if p^ = 2i/ then 


2z/- 1 

pP 


2z/- 1 

p/3 


= 1 > 


= 0 > 


'2v 

p/3 


'2v 

p/3 


- 1 . 


1 . 


Hence Vp{gcd{7i{Aj^i,),7i{Bj^i,)) < Vp{{2i/ — 1)!) which gives that 
gcd(7r(Hpj.),7r(Hj>)) | {2v - l)\ 

as required. 


□ 


It is worth remarking that for any fixed z/ > 4, we can End tighter 
lower bounds for the gcd by using Q directly. This can be used to 
tighten the results of Theorem HI for specific v. Unfortunately even 
when tightened in this way, we cannot achieve the conjectured bound. 
See Table H] 


We are now ready to prove Theorem 1.5 


Proof of Theorem 1.5. We fix integers u > 0 and 0 < j < k. We 
consider 2 cases: 
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If 0 < j < 1 / — 1 then we have 


e)fr)c:T'=e)(*r 


j J '{j + i^V- 


by our assumption on j weknow j + u < 2i/—1, hence (j + i^)! | (2z/—1)!, 
and therefore 

^k,j{k>) G (^qipyr^ ^ (2v-iy^- 
Otherwise we may assume that j > v. Then we have 

fk + v\ fk + v\ fk + 

{k + v) ■ ■ ■ {k + 1) ■ k ■ ■ ■ [k + V — i + 1) 


i! 

{k + v)---{k- i+ 1) 

z/! 

(i + 7! 

(A; + 7 ■ ■ 

■(A: + l) 

k ■ ■ ■ {k + V — j + 1) 

(A; + 7 ■ 

■ ■ (A: - i + 1) 


j! + 

(fc + z/) ■ ■ ■ (A; + 1) / /c 


(j + «^) ■ ■ ■ (j + 1) ■ J ■ ■ ■ (j - + 1) Vi - Vi 


z/! 



z/! 


Next observe that 

{k + p) ■ ■ ■ {k+ 1) f k V AV _ / k V/fc + z/V 

(i + «^) ■ ■ ■ (i +1) Vi - Vi/ Vi - Vi + 

(/c + z/) ■ ■ ■ (A; + 1) / A: \ A\ _ A + A\ 

+ V/-^/V/7 V / /V/7 

are both integers, hence there exists p, g G Z such that 


A 


k,j 


V = 


{k + v) ■ ■ ■ {k+ 1) 


(i + «^) ■ ■ ■ (i +1) ■ i ■ ■ ■ 7 - +1) \3 - ^ 
and where q | gcd((j + z^) ■ ■ ■ (i + I)) i ''' (i “ ^ + I))- 



I P I 

v\ = -v\ 

q 


2.2 


By Lemma 

there exists p' such that 


and transitivity of divisibility, q (2z/ — 1)! hence 


^fcjA) = 7 ■ 


(2Z.-1)!- 


Thus, for all integers z/ > 0 we have | ^ 2 v-i)\ desired. 


□ 
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3. A PROOF OF Theorem 


1.6 


(A‘) I r. 


Theorem 1.6 is an immediate corollary of: 


Lemma 3.1. Letp'^ > p°‘ and p'^ > u. Then the denom¬ 

inator of Apr_i ,^_i{y) is divisible by p'^. 

Proof. Let ■ Let p^ > p°‘ and p^ > v. Notice that 


Apr-l^y-liv) — 


P^ + V —1 

u-1 


p'^ - 1 

z/ - 1 


2 u-l 

v-l 


-1 


Consider the hrst term. 

''p^ -[- u — 1 
0-1 


(p" + n - 1) ■ ■ ■ (p" + 1) 


(z/-l)---l 

Observe that each factor of the top is eqnivalent mod p^ to the matching 
factor in the bottom. Hence ^ mod p. 

The second term is similar, with each term on the top eqnivalent 
mod p"^ to the additive inverse of the associated factor on the bottom. 
Hence = (—1)'^ mod p. 

Hence 

Apr-l,u-l{l^) - ^ 

with p co-prime to a. □ 

4. The case o = 3 and z/ = 4 

We see that ri = 1 = (j) and r 2 = 3 = ( 2 ). In this section we show 
the next two cases of Conjectnre 
and r 4 = 35 = (^). 

We hrst need the Lemma 


1.7 


hold, namely that ra = 10 = ( 3 ) 


Lemma 4.1. Let n and k be non-negative integers. If n is even and k 
is odd then Q) is even. 

Proof. By Knmmer’s theorem |H], 2 divides when there is at least 
one carry when k and n — k are added in base 2 . Since n is even and k is 
odd, n — k is, odd. The least signihcant bit of an odd integer represented 
in base 2 is always 1. Hence both k and n — k have a 1 in the least 
signihcant place. Thns when they are added, this will resnlt in a carry. 
So 2 divides (^). □ 

We now follow the proof of Theorem m nsing z/ = 3 to show: 


Theorem 4.2. Conjecture 1.1 holds for 0 = 3. That is r^ = ( 3 ) = 10. 
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Proof. We have that 10 |r 3 by Theorem 1.6 


As in the proof of Theorem 1.5, we hrst consider the case where 


0 < j < 2 . A quick calculation shows that 


Afc,o(3) 


10 

5{k + 3)k 
2 

{k — l){k + 2)(fc + 3)k 
4 


By considering the cases of k even or odd, we see that all of these values 
are always integers, and hence Akfl{3),Ak,i{3),Ak^2{3) G 
If j > 3 then, as in the proof of Theorem |1. 5, we have 


3 _ 3! fk + 3\ fk\ 

(J + W+2)(J + i)V J )\j) 

3! / k \ fk + 3\ 

^ j{j - l)(j - 2) [j - 3 ) Vj + 3/ 


We see that if 8 f gcd((j + 3)(j + 2)(j + 1), j{j - l){j - 2)) then 

2!3! 

Afcj(3) e 

as required. Hence we may assume that 8 | gcd((j + 3){j + 2)(j + 
1 ), j{j — l){j — 2)). Ifj is even then 8 | (j+ 3)(j+ 2)(j +1) implies that 
j = 6 mod 8 . We observe that 8 | j(j —l)(j—2) and 16 f j(j —l)(j —2). 
In this case we observe that one of (^-^ 3 ) and (j^^) is also even by 
Lemma 4T Hence we may write 

2 p 

AX3) = - A 

8 q 

where g is odd. This implies that 


2!3! 

Afcj(3) G 


as required. 

Similarly if j is odd, then j = 1 mod 8 , and 8 | (j + l)(j + 2)(j + 3) 
and 16 f (j + l)(j + 2){j + 3). Further one of and (^) is even, and 
hence 

^ , 2 p 
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where q is odd. Again this implies that 

2!3! 

Afcj(3) e 

as required. □ 


Theorem 4.3. Conjecture 1.7 holds for u = A. That is r^ = = 35. 


Proof. We have that 35|r4 by Theorem 1.6 


As in the proof of the previous theorem, we hrst consider the case 
where 0 < j < 3. A quick calculation shows that 


-4t,o(4)Q =36 

— 7k{k i) 

1 (ihb 7(k-l)k(k + 3){k + 4) 

- 

. (k-2)(k-l)k(k + 2)(k + 3)(k + 4) 

6lM(4)(^4j =-- 


By considering the various cases for k mod 12 (resp. 36), we see that 

these expressions are always integers, and hence Afc^o(4), Afc^i(4), G 

If J > 4 then, as in the previous proof, we have 


0+4)(i + 3)(i + 2)0 + l)V ] )[j) 

4! f k \fk + A\ 

j{j - l)(j - 2)(j - 3) \j - a) Vj + a) ■ 

From equation ([^ or Table we have that 
gcd((j + 4)(j + 3)(j + 2)(j + 1), - l)(j - 2)(j - 3)) I 7!/2 


Hence we have that Akj{A) G We still need to show that there 

is an additional factor of 3 in the numerator. 

To prove the result, we need to show that one of three things occurs 

• 9 tgcd((j + 4)(j + 3)(j + 2)(j + 1), j(j - l)(j - 2)(j - 3)) 

•3|(T)(T” 

If (j + 4) (j + 3) (j + 2) (j + 1) = j (j - 1) (j - 2) (j - 3) = 0 mod 9 then 
j = 2 mod 9 or j = 6 mod 9. Hence if j = 0,1, 3,4, 5, 7, 8 mod 9 
then Afcj(4) G as required. 
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k 

J 

a 

b 



k = Q mod 3 

j = 2 mod 3 

2 

2 

/ = 0 mod 3 

g = 1 mod 3 

k = 1 mod 3 

j = 2 mod 3 

4 

1 

/ = 0 mod 3 

g = 2 mod 3 

k = 2 mod 3 

j = 2 mod 3 

0 

3 

/ = 0 mod 3 

g = 2 mod 3 


Table 2. Cases when j = 2 mod 9 


k 

i 

a 

b 



k = Q mod 3 

j = 0 mod 3 

2 

4 

/ = 0 mod 3 

g = 2 mod 3 

k = 1 mod 3 

j = 0 mod 3 

1 

3 

/ = 0 mod 3 

g = 1 mod 3 

k = 2 mod 3 

j = 0 mod 3 

0 

2 

/ = 0 mod 3 

g = 2 mod 3 


Table 3. Cases when j = Q mod 9 


If j = 2 mod 9. then 27 f (j + l)(j + 2)(j + 3)(j + 4) so we have 
that 9 divides the gcd exactly. 

Consider 

fk + 4:\ fk\ _ fa,b{k,j) fk + a\ fk + b\ 

^ ^ [ J )[j) 9aAk:jA J )[ J J 

where fa,b{k,j) and ga,b{k,j) are polynomials. With careful choices of 
a and b we can construct fa,b and ga,b such that fa,b{k,j) will have more 
factors of 3 than gafi. 

For example, if a = 6 = 2 then 

f2,2{k,j) = {k + 4:){k + 3){k + 2 - j){k - J + 1) 

92,2ik,j) = {k - j + A){k - j + 3){k + 2){k + 1) 

Using the fact that j = 2 mod 3, we see that for k = 0 mod 3 that 
f 2 , 2 {k,j) = 0 mod 3 and g 2 , 2 {k,j) = 1 mod 3 and hence Q) = 0 
mod 3. A similar argument is given for k = 1 mod 3 and k = 2 
mod 3, summarized in Table Hence if j = 2 mod 9 then Afc j(4) G 
as required. 

If j = 6 mod 9 then 27 f j{j — l)(j — 2)(j — 3) so we have that 9 
divides the gcd exactly. 

Consider 

.ON ^ \ ^ la,b(J^^(k +a\fk + b\ 

Vi + 4/ \j - y 9a,b(k,j) Vi - 4/ Vi + V 

As before, we can break this into cases, as described in Table □ 


5. Additional support eor Conjecture 11.71 

We have computationally checked that for all k,j, v < 200 that Con¬ 
jecture 1.7 holds. Further, using the techniques of Theorems 4.2 and 
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4.3 we have computationally verified that for all j, z/ < 15 and all k that 


Conjecture 1.7 holds. It is not unreasonable to think that Conjecture 


1.7 can hold in general. Indeed, if we plot the non-integer entries in 


the lower triangular part of A{i') and colour them based on the prime 
factorization of their denominators in reduced form we obtain the frac¬ 
tal pattern seen in Figure ([^. This suggests that there is far more 
structure to the matrix A{i>) that we are currently exploiting. We note 
that from equation ([^ combined with Theorem 1.5 we would be able 
to prove that r^\2^ ■ 3^ ■ 7. We conjecture that = 2 ■ 3^ ■ 7. 

In this image of A(5), denominators are coloured red for 2, blue for 3, 
green for 7 and orange for 3^. If the denominator had contained any 
additional factors of 2, 3 or 5 then we would have coloured this value 
black. None occurred. Assuming that primes always give rise to the 
associated fractals early on, as seen in Figure [1} we would be led to 
believe that 4 f rs. 


6. Proof of Theorem 1.9: Wniv;2k) = n 


Proof of Theorem We rewrite (1.2) as 
Wn{m,2k) = 




k\ ■ {k + u)\ ■ z/!” ^ 




/ci! ■ ■ ■ kj. ■ {ki + iy)\--- {kn + p)! 


Let p = k he prime with 2z/ < p or let p = A; -|- z/ be prime with u < p. 
We claim that there does not exist indices 1 < i < j < n such that 
ki + v >p and kj + u > p. Indeed, this would lead to 

2p < ki + kj + 2v < (ki + ■ ■ ■ + kn) -|- 2z/ = /c -|- 2z/. 


\i p = k then 2z/ < p by assumption and hence 2p < fc -|- 2z/ < 2p, 
a contradiction. p = k + v then z/ < p by assumption and hence 
2p < (A; -|- z/) -|- z/ < 2p, a contradiction. 

If instead k = p^ and z/ = 0 it is easy to see that there does not exist 
indices 1 < i < j < n such that k^ + u > p^ and kj + u > p^. 

We consider 2 cases: 

If there exists 1 < i < n such that ki = k then clearly kj = 0 for 
j ^ i and hence 


A;! ■ (A: -|- z^)! ■ A"- ^ k\ ■ {k + u)] ■ AA ^ 

ki\ ■ ■ ■ kn ■ {ki + u)\ ■ ■ ■ {kn + z^)! A:! ■ 0! ■ ■ ■ 0! ■ (A: -I- z/)! ■ z/! ■ ■ ■ A 

Assume that ki < k for all 1 < i < n. 
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Figure 1. Non-integer entries of the first 1000 rows of A(5) 

If p = we see that p\k\ and p\{k -|- z/)!. We further see that at most 
one term in the denominator is divisible by p. Hence 

k\ ■ {k + z/)! ■ z/!”"^ 

/ci! ■ ■ ■ kj ■ {ki + v)\--- {kn + z/)! 

can be written as p| where p\ b, and thus is equivalent to 0 mod p. 

If p = k + u we see that p\{k + z/)!. We further see that no term in 
the denominator is divisible by p. Hence 

k\ ■ {k + z/)! ■ z/!”"^ 
ki\ - ■■ kj. ■ {ki + {kn + iy)\ 

can be written as p| where p\ b, and thus is equivalent to 0 mod p. 

If p^ = k and z/ = 0 we see that p'P^^\k\ and p'^^^\{k + v)\. We 
further see that we have at most 2p factors of p in the denominator, 
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with equality only if p\ki for all i. Hence 

k\ ■ {k + v)\ ■ 

fci! ■ ■ ■ kj ■ {ki + z/)! ■ ■ ■ {kn + z^)! 

can be written as where p f 6 , and thus is equivalent to 0 mod p^. 

Thus there are only n terms in the sum for Wn{i^', 2k) which are not 
0 mod p (resp 0 mod p^), namely when k^ = k for some k. In this 
case the term is 1 mod p (resp 1 mod p^) hence 

hTn(z/; 2k) = n mod p (resp. H 4 .( 0 ; 2k) = n mod p^) 

□ 


7. Comments 


We showed in Section that Conjecture |1.7| held for the case z/ = 
3 and z/ = 4. It is probably that this technique could be extended 
computationally for any hxed z/, although this is not clear. It is not 
clear that this technique would be extendable to arbitrary u without 
additional ideas. 

In Section we showed how the ideas of modularity of H7i(z/;/c) 
could be extended to fc = p^ or z/ > 0. It appears that something is 
also happening in the case when fc = p^ 7 ^ 4 and z^ = 1, although it is 
unclear how one would prove this. There are most likely many other 
relations that can be found when considering Wn modulo a well chosen 
prime power. 
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